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A convenient way to study phase transitions of finite spins systems of Unear size L is to fix 
boundary conditions that impose the presence of a system-size interface. In this paper, we study 
the statistical properties of such an interface in a disordered Potts ferromagnet in dimension d = 2 
within Migdal-Kadanoff real space renormalization. We first focus on the interface free-energy 
and energy to measure the singularities of the average and random contributions, as well as the 
corresponding histograms, both in the low-temperature phase and at criticality. We then consider 
the critical behavior of the interfacial adsorption of non-boundary states. Our main conclusion is 
that all singularities involve the correlation length (,av{T) ~ (Tc — T)~'^ appearing in the average 
free-energy F ~ {L / S,av{T))'^' of the interface of dimension ds = d — 1, except for the free-energy 
width AF ~ {L / S^variT))^ that involves the droplet exponent 6 and another correlation length 
£,var(T) which diverges more rapidly than S,av(T). We compare with the spin-glass transition in 
d = 3, where £,var{T) is the 'true' correlation length, and where the interface energy presents 
unconventional scaling with a chaos critical exponent C,c > ^jv [Nifle and Hilhorst, Phys. Rev. Lett. 
68, 2992 (1992)]. The common feature is that in both cases, the characteristic length scale Lch{T) 
associated with the chaotic nature of the low-temperature phase, diverges more slowly than the 
correlation length. 

I. INTRODUCTION 
A. Critical properties of interfaces in pure systems 

The critical points of statistical physics models are usually discussed in terms of bulk properties. However, it is 
also interesting to study how a critical system reacts when boundary conditions impose the presence of a system-size 
interface. For instance in an Ising ferromagnet defined on a cube of volume L'^, on may impose the (-I-) phase on the 
left boundary and the (— ) phase on the right boundary, with periodic boundary conditions in the other directions. 
The study of interfaces between coexisting phases near criticality has a long history [l| . The most important property 
is that the free-energy associated to the interface is proportional below Tc to its area with ds = d — 1 

Pure Ferromagnets : F"'*^''(L, T < T^) ^ h{T)L'^^ + ... (1) 

where the so-called interfacial tension fa{T) vanishes at criticality as a power- law [l| 

fo{T) ~ {Tc-rr (2) 

Finite-size scaling implies that near criticality, the interface free-energy F*"*'^'" should only depend on the ratio _L/^(T) 
between the system linear size L and the correlation length ^(T) oc (Tc — T)^" 



(i,T<T,)(x (^^y +... (3) 



This is similar to the requirement that the singular part of the bulk free-energy should scale as F**""^ {L/^{T)Y ~ 
(Tc - ry^L"^. The identification with the definition F''""'' - (T^ - Tf-°'L'^ in terms of the specific heat exponent a 
yields the hyperscaling relation 2 — a = vd. So the exponent ji of the interfacial tension fo{T) satisfies the Widom 
relation [ij 

li^vds = v{d - 1) (4) 

Exactly at criticality, the free-energy becomes of order punter f^j^^ ji^-j ^ 0(1), whereas the energy and entropy grow as 
E™*'^^{L, Tp) ~ L^^'^. Beyond these thermodynamic properties, the interest into interfaces was recently revived by the 
discovery |^] that some two dimensional critical interfaces are fractal curves which can be constructed via Stochastic 
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Loewner Evolutions (SLEs) reviewed in Q . Accordingly, the fractal dimensions of spin cluster boundaries of various 
two-dimensional spin models have been recently measured via Monte Carlo simulations in [1, Q • 

Whenever the system under study presents more than two phases, such as the Potts model considered in the this 
paper, a system-size interface between states 1 and 2 tends to produce a net adsorption of any non-boundary state, 
called state 3 here. This phenomenon of interfacial adsorption has been much studied in various pure models [6] with 
the following conclusions. The excess of state '3' due to the presence of a (1:2) interface with respect to the case (1:1) 
with no interface, defined as 

N„b = XI (< '^'^-3 >1:2 - < S„^^3 >l:l) (5) 
i 

is proportional to the area L''^ of the interface for T < Tc 

Nnt{L,T <Tc) ^woiT)L''' (6) 
Finite-size scaling argument yields that the coefficient wo{T) diverges at criticality as 

wo{T) cx in-Tf-' (7) 

where v is the correlation length introduced above, and where /3 is the order parameter exponent. This means that 
at criticality, the adsorption of non-boundary states Nnb scales as the global order parameter M — {Tc — T)^L'^. 

NntiL, Tc) ^ M{L, Tc) ^ L''-^ (8) 



B. Properties of interfaces below Tc in disordered systems 



In the field of disordered systems such as spin-glasses where the order parameter of the low-temperature phase is 
more complicated than in pure systems, it turns out that the properties of interfaces are very convenient to characterize 
the low-temperature phase via a so-called droplet exponent 'q 

Spin -glass: F'"*'=''(L, T < = T(T)L V + •- (9) 

where T(r) is a generalized 'stiffness' modulus and where up is a random variable of order 0{1). The exponent 6 is 
expected to satisfy the bound 6 < (d — l)/2 [7] The interface is expected to have a non-trivial fractal dimension dg 
with d — 1 < dg < d in the whole low-temperature phase [7] . This fractal dimension dg governs the energy and the 
entropy of the interface 

Spin -glass: E"''"'' {L,T < Tc) = a{T)L^ ue + ... (10) 

TS'^'^'^L, T <Tc) = a{T)L'^UE+ 

One actually expects the strict inequality 6* < ^ , so that the optimized free-energy of Eq [9] is a near cancellation of 
much larger energy and entropy contributions of Eq. [TlJ This is at the origin of the sensitivity of disordered systems 
to temperature changes or disorder changes, called 'chaos' in this context 0, M, H, E^j [U : roughly speaking, the 
chaos exponent ( = ^ — 6 governs the length scale L* ^ e^^/^ above which a small perturbation e in the temperature 
or in the disorder will change the state of the system. 

For non-frustrated disordered systems such as ferromagnetic spin models in dimension d, the interface below Tc is 
expected to be described by a directed manifold of dimension ds = {d — 1) in a random medium. In particular, in 
two-dimensional disordered ferromagnets, the one-dimensional interface is described by the directed polymer model 
[T^ . For this model, a droplet scaling theory has been developed [l3j | in direct correspondence with the droplet theory 
of spin-glasses [7| summarized above. In particular, the free-energy of the interface reads 

Random Ferromagnets : F™*"''(L, T < Tc) = fo{T)L'^^ + T{T)L^^'^^uf + ... (11) 

with a droplet exponent 9 which is exactly known to be 9{d — 2) ^ 1/3 on the two-dimensional lattice [T^. flsl. [T6| . 
The energy and the entropy of the interface reads [H, \vf\ 

Random Ferromagnets : £;™*'='^(L, T < = eo(T)L'*= + a{T)L^UE + ... (12) 

rS'™*"'^(L, T <Tc) = Tso{T)L'^'^ + a{T)L^UE + ... 

where the fluctuating term has again a bigger exponent ^ > 9 that the fluctuating term of the free-energy of Eq. [TT] 
As a consequence, the interface is again very sensitive to temperature or disorder changes with the chaos exponent 
C = ds/2 — 9. In particular in dimension d = 2, where the interface is a directed polymer in dimension 1 + 1, the 
chaos exponent is exactly known C = 1/2 — 1/3 = 1/6 [isl. ll§L fTol. [20j . 
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C. Properties of interfaces at criticality in disordered systems 

At criticality, the interface free-energy is expected to be a random variable up^ of order 0(1) 

F"*'='-(L,rc) + ... (13) 

For the spin-glass case, the interface free-energy of Eq. [9] is expected to scale as (L/^(T))^ in terms of the diverging 
correlation length ^(T) ~ {Tc — T)~'^, so that the critical exponent governing the vanishing of T(r) is Q 

T(T) - (Tc - T)"^ (14) 

which is the analog of Widom scaling relation for ferromagnets (Eq. |4|). 

For the energy of the interface, two possibilities have been described in the literature : 

(i) in the first scenario described in [7[ , the critical behavior follows the usual finite-size scaling forms in terms of the 
diverging correlation length ^(T) ^ {Tc — T)~'^. More precisely, the singular part of the energy or entropy is assumed 
to be of order l/(Tc — T) on the scale ^(T), so that the coefficient <j{T) in Eq. [TT] presents the following singularity 

1 / I \~ 

SG with 'Conventional' critical point : a(T) ~ — — I (15) 

Equivalcntly, one then obtains the following 'conventional random critical' behavior exactly at criticality Q 

SG with 'Conventional' critical point : E"'*'''{L, = L-ue, + - (16) 

where ue^ is a random variable of order 0(1). In our recent study of the directed polymer delocalization transitions 
on hierarchical lattices with 6 = 5 [2l| , we have found that the energy and entropy are governed by the 'conventional' 
critical behaviors of Eqs [T5l and UM 

(ii) however in (lol.lTl|. it has been found that a new exponent called the 'critical chaos exponent' can govern the 
response to disorder perturbations of spin-glasses at criticality, provided the inequality Cc > l/i' is satisfied. We refer 
to [lOl [TH for a detailed description of these chaos properties. Here, we will only mention an important consequence 
for the interface : it has been argued in [l^, [HI that this new exponent should govern the scaling of the interface 
energy at criticality 

SG with' Chaos critical exponent' : £;'"*^''(L, T^) = L'^'^ub^ + ... (17) 

in contrast with Eq. 1161 As a final remark on spin-glasses, let us mention that in = 2 where there is no spin- glas s 
phase (Tc = 0), recent studies have suggested that zero-temperature interfaces are actually described by SLE [2^l23||. 

For random ferromagnetic spin models, one expects 'conventional scaling' as in Eq. [16] for the energy where ue^ is 
a random variable of order 0(1). More generally, in the presence of relevant disorder, there is a lack of self-averaging 
in all singular contributions of thermodynamic observables in the sense that the leading term remains distributed 
[2i| . Note that for the Potts model with g > 3 states, the interface becomes a non-directed branching object at 
criticality. Some authors have studied the relevance of branching within a solid-on-solid approximation where the 
'directed' character of the low-temperature phase is kept [1^, H^. However, the 'directed' character is not expected 
to hold at criticality for at least two reasons : first, this 'directed' character does not hold at criticality already for 
pure ferromagnets, and second, in two dimensions, the directed polymer is always in its disordered dominated phase, 
whereas ferromagnets undergo a phase transition where the disorder relevance of the Harris criterion depends on q 
(see for a more detailed discussion). 

The aim of this paper is to study numerically the critical behavior of some two-dimensional random Potts ferro- 
magnet in the presence of relevant disorder. We have chosen to work on the diamond hierarchical lattice of effective 
dimension deff = 2, where large length scales can be studied via exact renormalization, and with the Potts model 
with q = 8 states so that disorder is relevant according to the Harris criterion (see Appendix |^ . We present detailed 
results on the statistics of the interface free-energy, energy, entropy and interfacial adsorption of non-boundary states. 



D. Organization of the paper 

The paper is organized as follows. In Section [TTl we recall the exact renormalization equations for the diamond 
hierarchical lattice, that are used to study numerically the disordered Potts model with q — 8 states on the diamond 
hierarchical lattice of effective dimension deff = 2. We then describe our numerical results on the interface free-energy 
statistics (Section IIII[) . on the interface energy and entropy statistics (Section IIV[) . and on the interfacial adsorption 
of non-boundary states (Section |V|. In Section IVTl we discuss the similarities and differences with the spin-glass 
transition in effective dimension de// = 3. Finally we give our conclusions in Section IVIII Appendix 1X1 contains a 
reminder on the pure Potts model on hierarchical lattices. 
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II. RENORMALIZATION EQUATIONS FOR SPIN MODELS ON HIERARCHICAL LATTICES 

A. Reminder on the diamond hierarchical lattices 




FIG. 1: Hierarchical construction of the diamond lattice of branching ratio b. 



Among real-space renormalization procedures [13] , Migdal-Kadanoff block renormalizations |28| play a special role 
because they can be considered in two ways, either as approximate renormalization procedures on hypercubic lattices, 
or as exact renormalization procedures on certain hierarchical lattices [l^, Is^] • One of the most studied hierarchical 
lattice is the diamond lattice which is constructed recursively from a single link called here generation n = (see 
Figured]): generation n ~ 1 consists of b branches, each branch containing 2 bonds in series ; generation n = 2 is 
obtained by applying the same transformation to each bond of the generation n = I. At generation n, the length L„ 
between the two extreme sites A and B is L„ = 2", and the total number of bonds is 

= (26)" = Li'"^'^ with deffib) = (18) 
where deff{b) represents some effective dimensionality. 



B. Spin models on hierarchical lattices 

On this diamond lattice, various disordered models have been studied, such as the diluted Ising model 'sT], ferro- 
magnetic random Potts model [13, [H, [s^l and spin-glasses [l3[ll[3i,[3i,[33,[3|. The random ferromagnetic Ising 
Hamiltonian reads 

Hjsing — ^ ^ ^ '^ij^i^j i^^) 
<ij> 

where the spins take the values Si = ±1 and where the couplings Jij are positive random variables (the spin-glass 
Hamiltonian corresponds to the case of couplings of random sign). The random ferromagnetic Potts Hamiltonian is 
a generalization where the variable ai can take q different values. 

H Potts = — ^ '2Ji,jScri,aj (20) 
<i,J> 

(We choose (2 J) to recover the Ising case for q = 2) 



C. Renormalization equation for the interface free-energy 



The free-energy cost F™^^"^ of creating an interface between the two end-points A and B of the diamond lattice of 
Fig. [T]is defined by 

pinter ^ - F++ (21) 
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where = —TlnZ^^ and = —TlnZ^ are the free-energies corresponding respectively to the same color at 

both ends or to two different colors at both ends. The renormalization equation are simpler to write in terms of the 
ratio of the two partitions functions and 



X = e 



z. 



++ 



z. 



and one obtains [32, [33, l34 



frf x'^'ht^ + iq-l) ' 



(22) 



(23) 



D. Renormalization equation for the interface energy 

The energy cost for creating an interface between the two ends at distance L is defined similarly by 

E-"" = -B,,^-^^ (1„ - in Z..) ^ I In |±± 
The renormalization equation reads in terms of the variable E = and x introduced above (Eq 

A x(^i)(a;(^^) - l)(x('^) +q- l)Et^ + x^'^\x^'''> - l){x^'^') +q- 1)eI^'^ 

En+l — 2^ 



1=1 



(a;(*i) + a;fe) + g _ 2)(x(^i)x(*2) + g - 1) 



(24) 



(25) 



E. Renormalization equations for the order parameter and the interfacial adsorption 

To study the order parameter and the interfacial adsorption, let us introduce the notation 

i 

for the number of spins in state '1' on a hierarchical lattice at generation n when the two end points A and B of Fig. 
[T]are respectively in states 'a' and 'b'. Using symmetries, one finally obtains closed renormalizations for the following 
five variables 



M. 



(1,1) _ 

n+l ~ 



M. 



(1,2) 
n+l 



E 

4=1 

b 

E 



i=i 



^n,n^n,.. (M^^^ + M^;^ -l) + {q- 1)(<-^) + M^;^) 



(27) 



(M^^ + M^:^ - 1) + a;„,, [M^J^ + M^^>) + {q- '^{M^^> + M^^>) 



(1,2) _^ , ,(2,2)s 



l-f(1.2) _L A/r(2,3)^ 



(2,1) 
n+l 



^j(2,2) 
-"'^n+l 



(2,3) 
n+l 



E 

i=i 

fc 

E 

1=1 

E 



(^4^^^ + - 1) + ^n,n=^n,.(i<-f + i<-f ) + (g - 2)(Mgf + 7i4^:f ) 

(A^S'^ + - 1) + (<i'^ + Ml^i^') + ^n,..(Ml':f) + M^f ) + (g - 3)(Ml^:f' + M^^f^) 



r(2,2 



g-2 

(2,1)n 



r,f(2,3) ^ , .(2a)^ 



The order parameter can then be defined as 

whereas the net absorption of non-boundary states of Eq. [5] reads 



(28) 
(29) 
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F. Numerical 'pool' method 



The numerical results presented below have been obtained with the so-called 'pool-method' which is very often used 
for disordered systems on hierarchical lattices : the idea is to represent the probability distribution P„ E„ ) of the 

interface free-energy F„ and energy at generation n, by a pool of N reahzations {{Fn^\En ^), .., {Fn'^\ En^^))}. 

The pool at generation {n + 1) is then obtained as follows : each new realization (i^^+i, £^^+1) is obtained by choosing 
{2b) realizations at random from the pool of generation n and by applying the renormalization equations given in Eq. 
[inland in Eq. [25l 

The initial distribution of couplings was chosen to be 

PpottsiJ) = 0{J>0)Je-'^ (30) 
for the ferromagnetic Potts case, and Gaussian for the spin-glass case 

Psg{J) = ^e-^ (31) 

At generation n = made of a single link (see Fig. [T]), the free-energy and the energy of the interface coincide and 
read in terms of the random coupling Ji drawn with either Eql^Ilor Eq. [3D] 

i^i2o = 4!lo - 2 J, (32) 

The numerical results presented below have been obtained with a pool of size N — 4.10^ which is iterated up to 
n = 60 or n = 80 generations. In the following sections, we study the random ferromagnetic Potts g = 8 on diamond 
lattice of effective dimension deff = 2 corresponding to a branching ratio b = 2. 



III. STATISTICS OF THE INTERFACE FREE-ENERGY 



As recalled in the introduction, the interface free-energy is expected to follow the scaling behavior of Eq[Tl1 below Tc 
and to become a random variable of order 0(1) at Tc (Eq |13p . In this section, we present numerical results concerning 
the singularities of the average and random contributions, as well as histograms, both in the low-temperature phase 
and at criticality. 

A. Flow of the average value and width of the interface free-energy 




In L Ln L 



FIG. 2: (Color online) Flows of the average value and of the width of the interface free-energy for many temperatures : (a) 
log-log plot of the average value F{L) of the free-energy distribution as a function of L. (b) log-log plot of the width AF{L) 
of the free-energy distribution as a function of L. 
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The flows of the average free-energy F{L) and of the free-energy width AF{L) are shown on Fig. [2] for many tem- 
peratures. One clearly sees on these log-log plots the two attractive fixed points separated by the critical temperature 
Tc- The value of Tc obtained by the pool method depends on the pool, i.e. on the discrete sampling with N values of 
the continuous probability distribution. It is expected to converge towards the thermodynamic critical temperature Tc 
only in the limit N — > oo. Nevertheless, for each given pool, the flow of free-energy allows a very precise determination 
of this pool-dependent critical temperature, for instance in the case considered 1.21685522 < T^""^ < 1.21685523. 

For T > Tc, both the average free-energy and the free-energy width decay exponentially in L. For T < Tc, the 
average free-energy grows asymptotically with the interface dimension ds = deff — 1 = 1 (see Eg [TT|) 

^(^) - (t-TF\^ + ••• '^itl^ = 1 (33) 



where ^av{T) is the correlation length that diverges as T — > T~ . The free-energy width grows asymptotically with 
the droplet exponent 9 (see Ea [TTt 



AF{L) ~ - — — with 0{b = 2) ~ 0.299 (34) 

\C,var\^ ) J 

where £,var{T) is the associated correlation length that diverges as T ^ T^T. Note that 9{b = 2) ~ 0.299 is the droplet 
exponent of the corresponding directed polymer model [H, [4^ . 



B. Divergence of the correlation lengths (,av{T) and ^var{T) 




The correlation lengths £,av{T) and (,var{T) as measured from the free-energy average value fEg [33)l and from the 
free-energy width (Eq [34] ) are shown on Fig. [3] (a). The log-log plot shown on Fig. [3] (b) indicates power-law 
divergences with two distinct correlation length exponents 

^av{T) ^^^{Tc-Ty'- with l^a. -1.07 

^var{T) ^oc^ {Tc ~ Ty with - 1-34 (35) 

In conclusion, our numerical results point towards the following singular behavior for the interface free-energy (see 
Eq. [11]) 

F"""iL,T < (^)" + {^)'u. + ... ,36) 
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where the average contribution and the random contribution involve two correlation lengths (T) and ^^ar (T) that 
diverge with distinct exponents fEq [55)1 . The presence of these two distinct correlation length exponents in the 
interface free-energy was a surprise for us, and we are not aware of any discussion of this possibility in the literature. 
The 'true' correlation length is expected to be ^av{T) that appears in the extensive non-random contribution to the 
interface free-energy. However, the presence of another length scale £,var{T) that diverges with a larger exponent 
remains to be better understood. 



C. Histogram of the interface free-energy below Tc 




-7 -6 -5 -4 -3 -2 -1 1 2 3 4 5 6 7 



FIG. 4: Statistics of the interface free-energy in the low-temperature phase (here T — 0.5) : Log representation of the 
distribution IIo{xf) of the rescaled variable xf = , which present asymmetric tails (see Eq. I38p 



In the low-temperature phase, the the interface free-energy is expected to follow the behavior of Eq 1111 where up 
is a random behavior of order 0(1). We show on Fig. |l]the probability distribution IIq{xp) of the rescaled variable 



XF 



F-F 
AF 



in log-scale to see the tails. The two tails exponents (77-, 77+) defined by 



lnno(a;F) ^, -\xf 

Xp — *-±oo 



(37) 



are compatible with the relations proposed in our previous work [sl] with ds — I, d — 2 and 9 ~ 0.299 (see Eg [34| 

d, 1 



V+ = 



1 



1.43 
2.85 



(38) 



D. Histogram of the interface free-energy at criticality 

At criticality, the interface free-energy is expected to become a random variable of order 0(1) (Eq. [T3)l : we show 
its probability distribution on Fig. [5] (a). To see the tails, we show in log-scale the distribution of the rescaled variable 
XF^^ on Fig. [51(b). 



IV. STATISTICS OF THE INTERFACE ENERGY 



In this section, we present the numerical results concerning the statistics of the interface energy. 
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FIG. 5: Histogram of the interface free-energy at criticality (a) Unrescaled probability distribution Pc{F) of the interface 
free-energy F at Tc (b) Log representation of the distribution I1c{xf) of the rescaled variable xp = 



A. Average and width of the interface energy 



As recalled in the introduction, the interface energy is expected to follow the scaling behavior of Eq [12] below 
Tc- The extensive non-random part eo(T) is directly related to the corresponding non-random part foiT) of the 
free-energy of Eq. [11] via the usual thermodynamic relation 

eo(r) = /o(T)-T^^ (39) 
As a consequence, the singularity found previously for fo{T) (Eq 

/o(T) ~ (-^) ° oc (T, - Tf^'-''- (40) 



U(T) 



determines the singularity of eo(T) near Tc 



eo(r) cx (Tc-T)^^"-^-' (41) 




FIG. 6: (Color online) Flow of the width AE{L) of the energy distribution as L grows (a) lnA_E(L) as a function of InL 
for many temperatures between T = 0.5 and T = Tc 4- e = 1.21685523 (b) Comparison of In A£;(L) (□), In AS(L) (0) and 
In A_F(L) (O) as a function of Ini at criticality. 
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We now consider the random contribution to the interface energy in Eq. [121 The flow of the width AE{L) as L 
grows is shown on Fig. [Hlfor many temperatures. For T < Tc, this width grows asymptoticaUy with the exponent 
ds/2 = 1/2 as expected (see Eq [T2|) 



AE{L) ~L^=L^ 

Exactly at criticahty, the energy width (and entropy width) grows as a power-law (see Fig[6]b) 

AE{L) ~ V-i" with ^ 0.92 
This value for j/c is in agreement with the value 1/vav — 0.93 (see Eq. l35| . 




30 





^fa, var 










(b) 





1.205 



1.215 



-5 -3 
In |T -T| 



(42) 
(43) 



FIG. 7: (Color online) Correlation lengths ^e{T) (*) and ^s{T) (0) as measured from the behavior of the energy and entropy 
widths (Eq. |45} (a) ln^i5(T) and \n^s{T) as a function of T, as compared to ln^a«(r) (□) and ln^.„ar(r) (Q) (b) In^siT) 
and ln^5(r) as a function of In \Tc — T\, as compared to In^au(r) and ln^i,ar(r) 



We now define a correlation length ^e{T) for T < Tc via the finite-size scaling form 

L 



AE{L,T) ~ L'^'^* 



(44) 



In the regime L ^ (^e{T), one should recover the L-dependence of the low-temperature phase of Eq[42l so the scaling 
function <i>(x) should present the asymptotic behavior $(a;) ~ 2;i/2-yc yielding the temperature dependence of the 
prefactor 



AE{L,T) 



i»fe(T) \ ^e{T) 



1 



l/2-ao 



L 



One similarly may define a correlation length ^s(r) from the finite-size scaling of the entropy width. 

As shown on Fig[7]b, the log- log plot presents some curvature, so that the asymptotic slope ve defined by 

^E{T)^{Tc-Ty^ 



(45) 



(46) 



is difficult to measure precisely. However, the slope i^e is close to the value i^av — 1-07 of Eq. [351 

In conclusion, our numerical results point towards the following singular behavior for the interface energy (see Eq. 

[El 



E'' 



\L,T<Tc) cc — — 
T-^Tr- Ic — ^ 



L 



£,av{T) 



ue 



(47) 



i.e. both the average contribution and the random contribution involve the same correlation length £,av{T). This 
result seems natural within the Fisher-Huse droplet theory 0, [l3| where the interface energy is a sum of random 
terms that follow some Central Limit asymptotic behavior. This picture is confirmed by the Gaussian distribution of 
the random variable ue that we now consider. 
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FIG. 8: Statistics of the interface energy in the low-temperature phase (here T = 0.5) : Log representation of the distribution 
Tio{xE) of the rescaled variable xe = ^^J? : it is a Gaussian 

B. Histogram of the interface energy below 

In the low-temperature phase, the the interface energy is expected to follow the behavior of Eq[T2l where ue is a 
random variable of order 0(1) which is expected to be Gaussian distributed within the droplet theory 0, [13 ■ this is 

in agreement with our numerical histogram of the rescaled variable xe = shown on Fig [51 

C. Histogram of the interface energy at criticality 




FIG. 9: Histogram of the interface energy at criticality (a) Unrescaled probability distribution Pn{E) of the interface free-energy 
E at Tc at generation n = 22 (b) Log representation of the distribution ric(a;_E) of the rescaled variable xe = ^^e 

We show on Fig. [3] our numerical results for the histogram of the interface energy at criticality : the unrescaled 
distribution for n = 22 generation is shown on Fig. [9] (a), whereas the distribution of the rescaled variable xe ^ ^KW 
is shown in log-scale on Fig. EKb). The fast decay of the tails show that the scaling of the interface energy at criticality 
is well measured via its variance (see Eq. I43p . 
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V. STATISTICS OF THE ORDER PARAMETER AND OF THE INTERFACIAL ADSORPTION 

A. Statistics of the order parameter M 

Below Tc, we find that the order parameter follows the scaling form 

M{L,T < Tc) ~ mQ{T)L'^ + mi(T)L'u, (48) 

The coefheient mo (T) of the extensive non-random term vanishes at criticality with the exponent 

moiT) (X (Tc - Tf with /3~ 0.163 (49) 

Exactly at criticality, the order parameter is expected to follow the behavior of Eq. [8l up to a random variable Uc of 
order 0(1) 

M{L,Tc)~L''-^Uc (50) 

Om previous measures of I'av — 1-07 (Eq. [35]) and /3 ~ 0.163 (Eq. [49]) would correspond with d = 2 to an exponent 
of order d — -J^ ~ 1.85. We measure (data not shown) 

M(L, T,) ^ AAf (L, Te) - L^-^^ (51) 
in agreement with the scaling relation of Eq ISOl 



B. Statistics of the interfacial adsorption Nnb 




(a) log- log plot of the average value N„b{L) as a function of L. (b) log- log plot of the width AN„b(L) as a function of L. 



The flows of the average value and of the width of the net absorption Nnb of non-boundary states (defined in Eq 
[5]) are shown on Fig. [TUl Below Tc, we find the scaling form 

NnbiT < Tc) - wq(T)L'^^ + 5{T)L^Vi with 4 = 1 (52) 

The coefficient wa{T) of the extensive non-random term is expected to diverge at criticality as in the pure case (Eq. 
[7]). The coefficient 5{T) of the random term is expected to diverge to yield the same finite-size scaling as Eq. [8]exactly 
at criticality, so that 

Nnb{L,Tc)^ L'-ivc (53) 

where Vc is a random variable of order 0(1). We measure (see Fig. [101) 

Nnb{L, Tc) - ANnb{L, Tc) ~ L^-^^ (54) 

again in agreement with the scaling relation of Eq [551 

The conclusion of this section is that our numerical results concerning the order parameter and the interfacial 
adsorption are consistent with 'conventional' scaling in terms of the correlation length ^ai,(T). 
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VI. COMPARISON WITH THE SPIN-GLASS ON DIAMOND LATTICE OF EFFECTIVE DIMENSION 

deff = 3 



The Migdal-KadanofF renormahzations with a branching ratio b = 4, which corresponds to an effective dimension 
deff = 3 (Eq[TH]), have been much used to study spin-glasses [13, UH, [H, [H, [13, [H] (for the case b = 2 corresponding to 
an effective dimension deff = 2 there is no spin-glass phase). As recalled in the introduction, a new 'chaos exponent' 

> l/i^ has been introduced in [To', Tl'l to describe chaos properties at criticality, and this exponent was argued to 
govern the energy at criticality fEq |17p . In the following, we confirm this scenario by directly measuring the statistical 
properties of the interface energy. We also discuss the similarities and differences with the random ferromagnetic case 
discussed above. 



A. Flow of the width of the interface free-energy 



15 



In 4F 



10 - 





10 20 30 40 50 60 

In L 



FIG. 11: (Color online) Spin-glass transition : log-log plot of the width AF{L) of the free-energy distribution as a function of 
L, for many temperatures. 



In the spin-glass case, there is no non-random leading term (Eq[H]) in contrast to the random ferromagnetic case ( 
Eg [TT|) . As a consequence, we only have to consider here the flow of the free-energy width AF{L) shown on Fig. [11] 
For T > Tc, the free-energy width decays exponentially in L. For T < Tc, the free-energy width grows asymptotically 
with the droplet exponent 9 (see Eq. [3]) 

/ r \ m 

AF{L) ~ — with 0{b = 4) ~ 0.255 (55) 

\C,varH ) / 

where £,var{T) is the correlation length that diverges as T — + T~ . The exponent 9{b = 4) ~ 0.255 is in agreement 
with previous measures [13, [H, [ll] . 

Again, the critical temperature obtained by this pool method depends on the pool, i.e. on the discrete sampling 
with N values of the continuous probability distribution. It is expected to converge towards the thermodynamic 
critical temperature Te only in the limit N oo. Nevertheless, for each given pool, the flow of free-energy width 
allows a very precise determination of this pool-dependent critical temperature, for instance in the case considered 
0.8810237 < TP°°' < 0.88102375. This value is in agreement with previous measures for a Gaussian initial condition 

The correlation length £,var{T) as measured from the free-energy width asymptotic behavior below Tc (Eq[55]) is 
shown on Fig. [12] (a). The plot in terms of the variable ln(Tc — T) shown on Fig. [12] (b) indicates a power-law 
divergence 

e™.(T)^a^(r, -T)-"^-- with i/™r^ 2.92 (56) 
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FIG. 12: (Color online) Spin-glass transition : Correlation lengths (,var{T) (O); ^e{T) (□) and ^s{T) (0) as measured from the 
behavior of the widths of the free-energy, energy and entropy distributions (a) ln^„ar(T'), ln^_B(r) and In^s(r) as a function of 
T (b) ln^i,ar(r), In^B(T) and ln^5(r) as a function of In \Tc — T\, : the asymptotic slope is the same and of order v^ar — 2.92 
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FIG. 13: (Color online) Spin-glass transition : Flow of the width AE{L) of the energy distribution as L grows (a) lnAE{L) 
as a function of Ini for many temperatures between T = 0.5 and T = Tc + e = 0.88102375 (b) Comparison of InAiJ(L) , 
In AS'(L) and In AF(L) as a function of InL at criticality {T^"^ = 0.8810237). 



B. Flow of the width of the interface energy 



The flow of the energy width AE{L) as L grows are shown on Fig. [13] for many temperatures. For T < Tc, this 
width grows asymptotically with the exponent ds/2 (see Eq. Ilip 

AE{L)~L^~L with = 4) = = 4) - 1 = 2 (57) 

Exactly at criticality, the energy width (and entropy width) grows as a power-law (see Fig [13] b) 

AE{L,Tc) ~ L'^" with Cc=i 0.58 (58) 

This exponent is clearly greater than l/uvar — 0.34 (see Eqs [56| and coincides with the chaos critical exponent 
measured in [Toj . 

We now define a correlation length £,e{T) for T <Tc via the finite-size scaling form 

AE{L,T) ~ L'^"^ ^ 

In the regime L ^ £,e{T), one should recover the L-dependence of the low-temperature phase of Eq[n71 so the scaling 
function $(a;) should present the asymptotic behavior <f>(a;) ^ a;^~^'= yielding the temperature dependence of the 
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prefactor 



AE{L,T) 



L 



As shown on Fig[T21 this leads to the same divergence as in Eg [551 

^E{T)^{T,-Ty--- with M„ar^2.95 



(60) 



(61) 



VII. SUMMARY AND CONCLUSIONS 



In this paper, we have studied the statistical properties of critical system-size interfaces in a disordered Potts 
ferromagnet. For the interface free-energy, our numerical results point towards the following singular behavior for the 
interface free-energy 



F™*"''(i,T < Tc 



fo{T)L' 



(X 



L 



T{T)UuF 

ds 



L 



^var (^) 



uf{T, L) 



(62) 



where the average contribution and the random contribution involve two correlation lengths S^av {T) and S^^ar (T) that 
diverge with distinct exponents at criticality fEq [35]) . The 'true' correlation length is expected to be £,av{T) that 
appears in the extensive non-random contribution to the interface free-energy. In particular, we have found that the 
interface energy follows the scaling form 



£;"'*-(L,T<r,) oc - 



1 



T 



L 



L 



Ue 



(63) 



i.e. both the average contribution and the random contribution involve the same correlation length (,av{T). This 
result seems natural within the Fisher-Huse droplet theory 0, [l3l | where the interface energy is a sum of random 
terms that follow some Central Limit asymptotic behavior. This picture is confirmed by the Gaussian distribution of 
the random variable ue, and by the 'conventional' behavior exactly at criticality 



L— 



UEa 



(64) 



However, the presence of another length scale S.variT) that diverges with a greater exponent Vyar > i^av remains to 
be better understood, in particular if one compares with the spin-glass transition. In the spin-glass case, S,variT) 
appearing in the random contribution of the free-energy is considered as the 'true' correlation length, since this is the 
leading term in the free-energy in this case. But then the interface energy is governed by some critical chaos exponent 
E'^"*^'^ {L,Tc) ~ L''" exactly at criticality, with C,c > ^/i^var, in contrast with the 'conventional' behavior of Eq. [M] So 
in both cases, even if the physical interpretation is different, one needs two different diverging length scales to describe 
the critical behaviors of the random contributions of the free-energy and energy or entropy. The physical origin seems 
to be in the chaos property of the random variable up{T, L) of order 0(1) in Eq. [SH Within one disordered sample, 
the random variable uf(T,L) strongly depends on the temperature, and this is why the energy and the entropy of 
the interface presents fluctuations that are not directly related to the scalings appearing in the free-energy. More 
precisely, the entropy can be obtained as a derivative of the free-energy with respect to temperature 



S'™*'='^(i,T < Tc) =- 



_ pi' 



\L,T) 



dT 
dh{T) 



dT 



dT{T) 



dT 



L%FiT,L) - T{T)L' 



, duF{T,L) 
dT 



(65) 
(66) 



(and similarly the energy reads E' — ^ ^ dT 

only of order L^, and thus we conclude that the fluctuation term of order a{T)L'^'^'^ with d 
entropy and in the energy ( Eq. I12p has for origin the derivative of the random variable uf{T, L). The identiflcation 
of these two terms yields 



The flrst term is the extensive term, the second term is 

d — 1 present in the 



T(r)L' 



,duF{T, L) 
df 



a{T)L^UE 



(67) 
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i.e. the derivative 



^^^W-T^r^) (69) 



T - Y(r) ^ UE^^^L<UE (68) 

is of order where = ^ — 9 is the chaos exponent of the zero-temperature fixed point. Let us now consider the 
singularity of the prefactor as T ^ Tc- If one defines a chaos length Lch{T) via 

one obtains that this chaos length Lch{T) diverges more slowly than the correlation length. More precisely, in the 
random ferromagnetic case, using T(r) - l/(^„ar(T))^ and a{T) - l/((Te - T)^iv''^ {T)), one obtains the singularity 

Random Ferromagnets : L,h{T) = Uy{T) x ( ^^^^ ) ' (70) 

So here the difference in scaling between the chaos length Lch{T) and the correlation length £,av{T) comes from the 
difference between ^yar{T) and ^av{T). In the spin-glass case, using T(r) ~ l/(^„ar(T))^ and cr(T) ~ {^var{T))'^''~'^^^'^ , 
one obtains the singularity 

Spin Glass : L,h{T) = Ur{T) x {Ur{T)) ^ (71) 

Here, the difference in scaling between the chaos length Lch{T) and the correlation length ^av{T) comes from the 
difference Q — > 0. 

In conclusion, beyond the differences in interpretation concerning the nature of the transition in disordered fer- 
romagnets (where the 'true' correlation length is ^av{T) associated to the non-random term of the free-energy, and 
where the interface energy presents conventional scaling at criticality E{L,Tc) ^ ii/'^-i-) and in spin-glasses (where 
the 'true' correlation length is Cvar(T) associated to the random term of the free-energy, and where the interface 
energy presents non-conventional scaling at criticality E{L,Tc) ^ L^" with a critical chaos exponent > ^/i^var), 
the common feature seems to be that in both cases, the characteristic length scale Lch{T) associated with the chaotic 
nature of the low-temperature phase, diverges more slowly than the correlation length. Note that for spin-glasses, 
Nifle and Hilhorst have found that the inequality C,c > l/i^ is satisfied in a finite range of dimensions d- < d < d+ 
above the lower critical dimension d^ , whereas for d > dj^, the usual scaling laws in terms of the correlation length 
exponent v are valid 10]. In disordered ferromagnets, one may similarly wonder whether the different singularities 
in (,av{T) and S,var{T) exist only in a finite range of dimensions d- < d < d^. It would be nice to clarify in which 
conditions the critical point of a disordered model is described by a single diverging length scale or by two diverging 
length scales. This probably requires a more precise understanding of the geometrical properties of the interface in 
the critical region that should be different for d < d+ and d > d^. 
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APPENDIX A: REMINDER ON THE PURE POTTS MODEL ON HIERARCHICAL LATTICES 

In the pure case, the renormalization of Eq. [23] reduces to the mapping T discussed in [s^, |4l| 



^"-^^(^"^n^trt^j =nx„) (Ai) 

The two attractive fixed points x = 1 (infinite temperature) and a; = 00 (zero temperature) are separated by a 
repulsive fixed point Xc (critical point). The critical exponents are obtained as follows [33l l4l| : the critical exponent 
v is determined by the linearized mapping around the critical point 
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and the specific heat exponent reads 



where def / — (ln(26) / ln(2) represents the dimension of the hierarchical lattice. In addition to usual power-laws, there 
are logarithmic oscillations coming from the discrete nature of the renormalization [4l| . For 6 = 2 (effective dimension 
deff =2), the critical point Xc(q) corresponds to 

q = {x,-l){V^,-l) (A4) 

and the transition is second order for any q. The Harris criterion indicates that disorder is relevant for 

q > qHarr^s{b = 2) = 4 + 2^2 - 6.828... (A5) 

and this is why we have chosen to use the value q = 8 for our numerical simulations of the disordered case. 
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